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Abstract. Differential rotation is known to suppress linear instabilities in fusion 
plasmas. However, numerical experiments show that even in the absence of growing 
eigenmodes, subcritical fluctuations that grow transiently can lead to sustained 
turbulence, limiting the ability of the velocity shear to suppress anomalous transport. 
Here transient growth of electrostatic fluctuations driven by the parallel velocity 
gradient (PVG) and the ion temperature gradient (ITG) in the presence of a 
perpendicular (E x B) velocity shear is considered. The maximally simplified (but 
most promising for transport reduction) case of zero magnetic shear is treated in 
the framework of a local shearing box approximation. In this case there are no 
linearly growing eigenmodes, so all excitations are transient. In the PVG-dominated 
regime, the maximum amplification factor is found to be with N cx q/e (safety 
factor/aspect ratio), the maximally amplified wavenumbers perpendicular and parallel 
to the magnetic field are related by kyPi « (e/q)^/^fc||Wthi/5', where pi is the ion 
Larmor radius, Vthi the ion thermal speed and S the E x B shear. In the ITG- 
dominated regime, N is independent of wavenumber and N cx Vt\iil {LtS), where Lt 
is the ion-temperature scale length. Intermediate ITG-PVG regimes are also analysed 
and N is calculated as a function of q/e, Lt and S. Analytical results are corroborated 
and supplemented by linear gyrokinetic numerical tests. Regimes with N < 1 for all 
wavenumbers are possible for sufficiently low values of q/e (< 7 in our model); ion-scale 
turbulence is expected to be fully suppressed in such regimes. For cases when it is not 
suppressed, an elementary heuristic theory of subcritical PVG turbulence leading to 
a scaling of the associated ion heat flux with q, e, S and Lt is proposed; it is argued 
that the transport is much less "stiff" than in the ITG regime. 
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1. Introduction 

It has long been understood that anomalous transport in tokamaks is caused by 
turbulence at or just above the ion Larmor scale, which is powered by drift instabilities 
extracting free energy from background gradients. The ion-temperature-gradient (ITG) 
instability [lOl |15l [31] and the resulting turbulence [H] have been identified as a 
particular culprit. With the advent first of gyrofiuid and then gyrokinetic numerical 
simulations, ITG turbulence has been the main focus of numerical studies aiming 
to map out levels of turbulent transport expected for any given set of equilibrium 
parameters (mean profile gradients and the magnetic configuration) [30l [181 EHl [H]- 
These studies assumed that the differential (toroidal) rotation of the tokamak plasmas, 
caused by momentum injection from the heating beams, had a negligible effect on 
the turbulence and on the resulting transport. This, however, is known to be an 
inadequate approximation for many devices and configurations in which the rotational 
shear can be comparable to the typical turbulent rate of strain. Linear eigenmode 
analysis [3l [5l [201 E] suggests that the ITG growth rates are reduced and can be 
completely quenched hy E ^ B shear. A reduction or even complete suppression of 
associated transport was therefore expected and indeed found numerically, at least for 
some parameter regimes [281 [33 IS] — a result that kindled high hopes for controlling 
turbulence with shear and achieving transport bifurcations to steeper temperature 
gradients [211 HH] • 

Since the tokamak rotation is (to lowest order in the gyrokinetic expansion) purely 
toroidal [271 HSl HOI UHl [2], strong perpendicular E x B shear comes at the price of 
a (stronger by a factor of q/e) parallel velocity gradient (PVG) — which is itself a 
source of free energy and so can trigger a drift instability which in turns gives 
rise to turbulence and anomalous transport [l?! [191 U] • This instability is in fact also 
suppressed by the E x B shear, but only in the sense that no unstable eigenmodes 
survive at large enough velocity and small enough magnetic shear: transient linear 
amplification is still possible [45l \M\ [35]. This transient amplification (illustrated in 
figure [1]) can be sufficient to give rise to subcritically excited turbulence [3 [251 [26] . 
It is this turbulence that limits the effectiveness of differential rotation in suppressing 
anomalous transport. While transport bifurcations are still possible [251 [261 [37] , finding 
them in the parameter space turns out to be quite a delicate task, mostly because the 
subcritical PVG turbulence is a largely unexplored phenomenon. In particular, it is not 
amenable to the usual "quasilinear" mixing-length-type arguments because those are 
based on the calculation of linear growth rates — and it is not obvious what should 
replace them for subcritical turbulence. 

We take the rather obvious view that addressing this problem should start 
from a systematic understanding of the linear transient amplification of gyrokinetic 
fluctuations in the presence of velocity shear. The purest example is presented by 
the limit of zero (negligibly low) magnetic shear, which has the twin advantages of 
analytical simplicity and of being most amenable to transport bifurcations or at least 



Suhcritical fluctuations in rotating plasmas 3 
le+04 



le+03 

X 

^ le+02 
I le+01 
le+00 

100 200 300 400 500 600 

t/ I Vthi 

Figure 1. Time evolution of the squared amplitude (normalised to its initial value) of a 
pure-PVG-driven linear perturbation, obtained in a direct linear numerical simulation 
using the gyrokinetic code AstroGK [36] with 1/Ln = 1/Lt = 0, q/e = 50, t/Z = 1, 
kyPi — 1, k\\Vthi/S — 0.5. The effective growth rate for this case is shown in figureH) 
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reduced levels of turbulent transport — as suggested both by numerical experiments 
[25| |26] and laboratory measurements [32l |33]. Formally, this regime supports no 
growing eigenmodes at any finite value of velocity shear, so we can focus on transient 
amplification and its dependence on the parameters of the problem (the ion temperature 
gradient, the velocity shear, the safety factor q and the device aspect ratio e) without 
the complications of dealing with the transition to linear eigenmode stability. This is 
the kinetic treatment of this problem, following in the footsteps of the fluid theory [35] 
(which was done for a finite magnetic shear; note that zero magnetic shear is formally 
a singular limit). 

Since much of the transport modeling for fusion devices relies on gyrokinetic 
simulations, it is important to ascertain that the linear behaviour predicted analytically 
is reproducible in direct gyrokinetic simulations with standard numerical codes. In 
what follows, we do this using the AstroGK code [36], which, for all practical purposes, 
is the slab version of the widely used fusion code GS^ (running GS2 itself in slab mode 
produces similar results). The code solves (j2]) and ([3]) with nonlinearity switched off 
and a collision operator jll |6] added to provide small-scale regularisation in phase space. 
These simulations demonstrate the degree to which our asymptotic results represent an 
adequate description of realistic parameter regimes and how the transition from the 
short-time to the long-time limit occurs, subject also to (minor) modifications of our 
results by collisions, finite resolution and numerical inaccuracies. 

The rest of the paper is organised as follows. In §[2|, we introduce gyrokinetics in 
a shearing box — the governing equations of the minimal model we have chosen to 
treat. In §[3l we consider a further simplified limit in which the ITG drive is negligible 
compared with the PVG drive and so we can concentrate on the essential properties of 
the latter, namely, work out how the transient amplification time and the amplification 
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exponent depend on the velocity shear, q and e and which wavenumbers prove most 
prone to being amphfied. In §|H we generahse these considerations by including the ITG. 
Both analytical and (linear) numerical results are presented throughout. A qualitative 
summary of the linear regime and a comparison with the fluid limit treated in [35] are 
given in §[5l Finally, in §[6], a criterion for the onset of subcrtitical PVG turbulence is 
proposed, followed by a very crude heuristic theory of this turbulence and of the resulting 
heat transport — these considerations provide a version of the standard mixing-length 
"quasilinear" arguments suitable for transiently growing fluctuations and a set of scaling 
predictions in the spirit of [8]. In particular, we give a semi- quantitative form to the 
argument that heat transport must be much much less "stiff" in the presence of velocity 
shear [321 |33l [251 126] than in in the standard ITG regime because the turbulence is no 
longer driven by the temperature gradient. 



2. Gyrokinetics in a shearing box 

Consider a plasma in a strong mean magnetic field. On the assumption of axisymmetry, 
this field can be expressed in the form B = x V0, where (magnetic flux) 

and are two scalar functions determined by solving the mean MHD equilibrium 
equations and is the azimuthal angle with respect to which symmetry is assumed. It 
can be shown (see [2] and references therein) that if such a plasma rotates with some 
velocity ordered in the gyrokinetic expansion as similar in size to the speed of sound, 
this mean rotation must be purely azimuthal, the same for both species and with an 
angular velocity co that depends on the flux label ip only: u = uj{tp)R'^'V(j), where R is 
the radial coordinate. Thus, each of the nested toroidal flux surfaces rotates at its own 
rate and there is a velocity shear set up by the variation of u with ip. In this paper, 
we will be concerned only with the effect on the plasma stability of this differential 
character of the rotation. Formally, this effect is isolated by assuming the Mach number 
M = u/vthi {vthi is the ion thermal speed) to be moderately low. Mathematically, this 
can be cast as an expansion in M <^ 1 (subsidiary to the gyrokinetic expansion) where, 
while the mean velocity is ordered subsonic, u ~ 0{M), its scale length is ordered as 
0{1/M). Under this scheme, effects such the Coriolis or centrifugal motion (as well 
as other, more obscure, ones) that scale with the magnitude of the mean velocity are 
negligible, while velocity gradients are retained. 

Let us consider the vicinity of some flux surface ip = ipo and introduce a local 
orthogonal Cartesian frame (x, y, z) that is moving with the flux surface ipQ and in which 
X is the cross-flux-surface ("radial") coordinate, z is the coordinate in the direction of 
the magnetic field and y completes the orthogonal frame [y = z^x). Then the velocity 
field can be replaced by a pure linear shear flow: 

«.V«S.^. S=^-, (1) 



where Bp = \'Vijj\/R is the poloidal magnetic field (see Appendix A.l for details; figure 



|2] illustrates the geometry). This model might be called the "flying slab approximation," 
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or, perhaps more conventionally, the shearing box — a common analytical simplification 
in the fiuid dynamics of differentially rotating systems [23] . 

We will make three further simplifying approximations by assuming all fluctuations 
to be electrostatic, electrons to have purely Boltzmann response and the mean magnetic 
field to be (locally) straight and uniform, so it has neither curvature nor variation of 
magnitude nor shear (which means that the magnetic drifts can be dropped). With 
all these assumptions and working in the fiying slab (shearing box), we can write the 
gyrokinetic system of equations in the following form (see Appendix A ) 

/ d{ip)R dh d{ip)Rdh^ 
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where h{t, R, w, w\i) is the ion gyrocentre distribution function, Fo^ip, w) the equilibrium 
Maxwellian (of the ions), (p{t,r) = Zecp/Ti the nondimensionalised electrostatic 
potential {Ze is the ion charge), pi the ion Larmor radius, rii the ion number density, w 
the particles' peculiar velocity with respect to the mean fiow, r the position coordinate, 
R the guiding-centre coordinate, and the angle brackets denote the gyroaverages holding 
constant the coordinate appearing in their subscript (the precise definitions are provided 
in [Appendix A[). We have defined the equilibrium gradients and two other standard 



parameters as follows 
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(4) 



where Ui and Tj are the ion density and temperature, respectively, and i?^ is the 
azimuthal mean magnetic field. 
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Note that the velocity shear S appears in two ways in as perpendicular shear 
(multiplying xd/dy) and as parallel shear (multiplying Fq in the right-hand side). The 
perpendicular shear rips apart the unstable fluctuations and will have a stabilising 
effect, whereas the parallel shear (the PVG) acts as a source of free energy in a manner 
analogous to ITG and drives a drift instability (§ l3.ip . The relative size of these two 
effects is set by the value of q/e. 

2.1. Case of non-zero magnetic shear 

Including a (locally) constant linear magnetic shear into the problem amounts to 
replacing in ([2]) 

dh ( ^ w\\\ dh 



where Lg is the scale length associated with the magnetic shear. This appears to 
introduce complications as we now have an "effective shear" that depends on the particle 
velocity w\\ . However, since the size of w\\ is constrained by the Maxwellian equilibrium 
distribution, this term can be neglected provided Mg = SLs/vtu ^ 1. Under this 
assumption, the theory developed below applies without modification. We note that 
the "shear Mach number" is precisely the parameter that is known to control linear 
stability and transient amplification of in the ITG-PVG-driven plasmas in the fluid 
(collisional) limit [35] . 



2.2. Shearing frame 

The next step — standard in treatments of systems with linear shear — is to make a 
variable transformation (t, r) — {t',r') that removes the shear terms {Sxd/dy): 

t' = t, x' = X, y' = y — Sxt, z' = z, (6) 

and similarly for (t, R) — ^ {f, R'). The Fourier transform can then be performed in the 
primed variables, so 

V = J2^k'{ty'^'-^' = J]^M'(^')e^'^'''^■^ (7) 
fc' fc' 

where = k'^ — Sk'yt', ky = k'^, and /cy = fcy (we denote = k^). As usual in the 
gyrokinetic theory, working in Fourier space allows us to compute the gyroaverages in 
terms of Bessel functions: 

= E ^o(a(t'))^fc'(i')e^''-''', {h)r = Yl ^o(a(t'))/^fc'(i')e^''"', (8) 
fc' fc' 

where a{t') = k^w^/Vti = {w^/Vti)^ {k'^ - Sk'^Vf + k'!^ = {k'yW±/ni)Vl + SH"^, and 
we have shifted the origin of time: t" = t' — S~^k'^/ky. 

Finally, we rewrite the gyrokinetic system (I2])-(I3]) in the new variables {t", k'). Since 
we are interested only in the linear problem here, we will drop the nonlinearity. We also 
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suppress all primes in the variables. The result is 

3 
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FM<t))^u, (9) 
(10) 



where = kyPiVthif^Ln is the drift frequency and rji = Ln/ Lt] the argument of the 
Bessel function is a{t) = {w±/vthi)kyPi\/l + SH^. 

2.3. Integral equation for the linearised problem 

We integrate ^ with respect to time, assume the initial fluctuation amplitude small 
compared to values to which it will grow during the subsequent time evolution, rescale 
time \k\\ |fth«^ t, denote At = t — t' , use fllOl) . in which the velocity integrals involving 
the Maxwellian Fq = ne~"'^''^thi/(7ri;^j^j3/2 done in the usual way, and obtain finally 
the following integral equation for <^(t): 



(l + ^ - ro(A, A)) m = I dAt e-^*^/^ { (' - l) f 



— «CJ* 



1 + A(A,A0-1 
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where cj* = uj*/\k\\\vthi = kypi/2\k\\\Ln is the normalised drift frequency, us 
Skypi/k\\Vt\ii the normalised shear parameter, and 



ro(A, A') = e-(^+^')/2/o(v^), A(A, A') 



^_A + A:^^/.(v/AA^) 



(12) 



2 /o(VAAO 

where \{t) = {klpf + u)lt^)/2, A' = A(t') = (fc^pf + ult'^)/2, and Jq and h are modified 
Bessel functions of the first kind. 

Equation (|TT|) is the master equation for the linear time evolution of the plasma 
fluctuations driven by the ITG (the w,,,?], term) and the PVG (the {q/e)uJs term). 



3. Solution for the case of strong shear 

We will first consider the maximally simplified case of pure PVG drive (strong shear). 
This is a good quantitative approximation to the general case if a)*, r^jo;* ^ {q/e)ujs, 
which in terms of the basic dimensional parameters of the problem translates into 

qS vthi vthi 

This is the regime into which the plasma is pushed as the flow shear is increased — 
under certain conditions, the transition can occur abruptly, via a transport bifurcation 
[25| l37t 126] . Besides being, therefore, physically the most interesting, this limit also has 
the advantage of particular analytical transparency (the more general case including 
ITG will be considered in §|1]). 



Suhcritical fluctuations in rotating plasmas 



8 



Thus, neglecting all terms that contain and ?7j, fllip becomes 
(l + ^ - ro(A, A)) m = dAtAte-^''/' {^^cos - l) ro(A, X'Mt - At). (14) 

3.1. Short-time limit: the PVG instability 

Let us first consider the case in which the velocity shear is unimportant except for the 
PVG drive, i.e., we can approximate A ~ A' ~ kyp1/2 and so there is no time dependence 
in the Bessel functions in ( !T4|) . We would also like to be able to assume t S> 1 so that the 
time integration in f ll4p can be extended to oo. Formally this limiting case is achieved 
by ordering kyPi ~ 1 and l^t^ 1/us ^ q/e (in dimensional terms, this is equivalent 
to St ^ 1, k\\Vthit ^ 1 and 1 <^ k\iVthi/S < q/e). Physically, this regime is reahsed in 
the initial stage of evolution of the fluctuations or, equivalently, in the case of very weak 
shear but large q/e. 

Under this ordering, we can seek solutions to (fn|) in the form (f(t) = (poexp^—iut), 
where u = |fthi is the nondimensionalised complex frequency and u its dimensional 
counterpart. The time integral in 014p can be expressed in terms of the plasma dispersion 
function, which satisfies 



poo 

Z{cu) = I \ rfA^e''^^*-^*'/^ Z\uj) = -2[l+cuZ{uj)]. (15) 
Jo 

With the aid of these formulae, (ITT]) is readily converted into a transcendental equation 
for Cj: 

^^^'^ ' ^^a;5-l)[l + a;Z(a;)], (16) 



ro(A) ve 

where ro(A) = e~^/o(A), A = and Gjs = SkyPi/k\\Vtu- 

Equation ( fT6i) is simply the dispersion relation for the ion acoustic wave modified 
by the PVG drive term. This point is probably best illustrated by considering the cold- 
ion/long- wavelength limit r < 1, u; = u;/|A;|| |i;thi > 1, A = kyp^ /2 ^ 1. Then ro(A) ^ 1, 
1 + QZ{u)) ~ — l/2a;^ + iGj-Jii e^^^ and so, restoring dimensions in f lT6|) . we get 

1/2 



f (l-^c..) ^ ..±.„c.(l-^i^) , (17) 

where in the last expression, we have restored dimensions and denoted = 
(Z/2r)"'^/^fthi = {Te/TniYf^ , the sound speed. When {q/e)oJs is sufficiently large, the 
sound wave is destabilised and turns into the PVG instability. Note that it loses its real 
frequency in this transition. 

Some further (elementary) analytical considerations of the PVG instability are 



relegated to Appendix B.l Here it will suffice to notice that if the (dimensional) growth 
rate 7 = Imcj is scaled by qS/e and /cy by {qS/e)/vthi, their mutual dependence is 
universal for all values of the velocity shear or q/e, namely, 

7 — /I f^yPi 



e \ qS/e 




Figure 3. Left panel: Growth rate 7 = Imtj (normalised by qS/e) vs. k\\Vthi/{qS/e). 
Here t/Z ~ 1, and the two curves are for kyPi = 0.1 (red) and kyPi = 1 (brown). 
The growtli rate becomes independent of ky for kyPi ^ 1 (see end of [Appendix B.l[ ) 
and all curves for large kyPi are very close to each other and similar in shape to 
kyPi = 1 (see right panel). The mode has no real frequency for k\\vthi/iqS/e) < kyPi] 
for k\\Vt\^i / {qS / 1) > kyPi, it turns into a damped sound wave: the corresponding 
frequencies and the damping rate are shown only for kyPi = 0.1, as thin blue (dashed) 
and red (solid) lines, respectively. The discrete points show growth rates calculated 
by direct linear numerical simulation using the gyrokinetic code AstroGK |36j . Right 
panel: Contour plot of "//(qS/e) vs. k^\Vthi/{qS/e) and kyPi. Only positive values are 
plotted, black means 7 < 0. The red curve shows the stability boundary (jB.ll) . The 
white line shows the boundary (jB.2[) between PVG modes (above it) and sound waves 
(below it). 



Once this rescaling is done, the dispersion relation f lT6|) no longer contains any 
parameters (except r/Z, which we can safely take to be order unity). The growth 
rate, obtained via numerical solution of ( fT6l) with r/Z = 1, is plotted in figure El The 
maximum growth rate is 7max ~ 0.10(gS'/e). This peak value is reached when kyPi ^ 1.0 
and k\\ ^ 0.10(gS'/e)/fthj, i.e., at {q/e)us ~ 10. 

The conclusion is that, at least in the initial stage of their evolution, plasma 
fluctuations in a significant part of the wavenumber space {ky and k\\ are constrained 
by ( IB.ip ) are amplified by the PVG. This amplification does not, however, go on for a 
long time as the approximation we adopted to derive the dispersion relation ( !T6|) breaks 
down when ust ~ kyPi (or St ~ 1 if the dimensional units of time are restored) — this 
gives 7t ~ 0.1(g/e), so, realistically, after barely one exponentiation. The key question 
is what happens after that. We will see shortly that all modes will eventually decay and 
that the fastest initially growing modes are in fact not quite the ones that will grow the 
longest or get maximally amplified. 
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3.2. Long-time limit: transient growth 

Let us now investigate the long-time limit, w^t ^ 1, kyPi (or, in dimensional form, 
St ^ 1, kyPiSt ^1). In this limit, the kernel involving the Bessel function in ( |T4l) 
simplifies considerably: we have A ~ Ugt'^/2 ^ 1, A' ~ tij|(t — At)^/2 ^ 1 and so 

-{v^-v^)2/2 -4Ai2/4 

Working to the lowest nontrivial order in 1/t, we can now rewrite f|T^ as follows 
(l + ^) \us\tm = ^ l^^dAtAte-^^'i+mty^ (^l^^ _ _ A^). (20) 

We will seek a solution to this equation in the form 



ip{t) = v^oexp 



dt'^{t') 



(21) 



where 7(t) = 7(t)/|fc|||fthi is the effective time- dependent growth rate (nondimension- 
alised) and ■jit) the dimensional version of it (remember that time is scaled by |A;|||fthj)- 
Because of the exponential in the kernel, the memory of the time-history integral in the 
right-hand side of (!20|) is limited, so At <^ t and we will be able to make progress by 
expanding 

"* At^ 
dt'j{t')-Atj{t) + —j'{t) + .. 



(p{t — At) = (fo exp 







(22) 



We will assume that this expansion can be truncated; the resulting solution will indeed 
turn out to satisfy 7'(t) <^ 1, with all higher-order terms even smaller. Substituting 
fl22|) into fl20|) . we obtain an implicit transcendental equation for 7(t): 

(l + L)lajs\t = ^ r dAt At e-A*7(t)-{4+i)AtV4 f ^ - 1 V (23) 

This equation can be written in a compact form by invoking once again the plasma 
dispersion function f|T5|) : denoting jit) = j{t)/ a/u;| + 1, we get 

(l + ^) ^{ul + l)\us\t =[-^0Js- l) [1 + il{t)Z{iim (24) 

— effectively, a time- dependent dispersion relation, reminiscent of the PVG dispersion 
relation (fT6|) . 

3.2.1. Transient growth. First of all, it is immediately clear from that as time 
increases, (the real part of) the effective time-dependent growth rate j{t) must decrease 
and indeed go negative because the right-hand side has to keep up with the increasing 
left-hand side. Therefore, fluctuations will eventually decay. However, if Re7(t) is 
positive for some significant initial period of time, there can be a substantial transient 
amplification. 

We can determine the time to when the transient growth ends by setting Re 7(^0) = 
in We immediately find that Im7(to) = as well and that 

{q/e)Qs - 1 .2,^ 
° {l + r/Z)^{ujl + l)\us\' ^ ^ 
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Figure 4. Time evolution of the effective growth rate: the red (bold) line is 
7(t) = 7(t)/|fc|| |uthi\/w| + 1 obtained as a numerical solution of ((26)) and plotted vs. 
t/to (the time axis is logarithmic in base 10). The black (thin) lines are the asymptotics 
([27]) and (|32|) . The discrete points show the time evolution of the effective growth rate 
obtained in a direct linear numerical simulation using the gyrokinetic code AstroGK 
[36] with 1/L„ = 1/Lt = 0, q/e = 50, t/Z = 1, kyp, = 1, k\\VtM/S = 0.5. The 
short-time-limit PVG growth rate for this case (obtained by solving (|16|) ') is shown as 
a dotted horizontal line. The time is normalised using ([25]) for t^. The time evolution 
of the perturbation amplitude for this case is shown in figure [TJ 



The dependence of to on ky and k\\ — via us = SkyPi/k^i^Vt^i and via the time 
normahsation factor of |A;||fthj| — tells us which modes grow longest. The interesting 
question, however, is rather which modes get maximally amplified during this transient 
growth. 

3.2.2. Maximal amplification. The total amplification factor is given by e^, where 
= dt'y{t) is the number of exponentiations experienced by the mode during its 
growth period. In order to determine this, we need to know the time evolution of 7(t) 
up to t = to- Using ( l25l) . it is convenient to rewrite (HM as follows 

- = l + i^Z{i^). (26) 

^0 

When t to, the solution is found by expanding the plasma dispersion function in 
7 > 1: 

t I ft^ 

This asymptotic is not valid when t approaches to- More generally, ( 126|) has a solution 
7 = 7(t/to), whose functional form is independent of any parameters of the problem. 
It is plotted in figure S] together with the asymptotic (1271) and with a direct numerical 
solution showing how the transition between the short-time (§ 13. ip and long-time limits 
occurs. The amplification exponent is easily found: 

N = r dt^it) = toM+~l f ^ 0.45 7—^^^^%r^j--. (28) 

Jo ^ Jo {1 + T/Zjy/Loi; + l\UJs\ 
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Figure 5. Left panel: The amplification exponent N vs. k± and k\\ in tlie same 
AstroGK simulation as used to produce the discrete points in figure 2] The straight 
line shows the relationship between the wavenumbers given by the second formula 
in (j29l) . Rigtit panel: The maximum amplification exponent A^max [N given by (j28p . 
maximised with respect to Cos] vs. q/e. The dotted lines show the g/e 3> 1 asymptotic 
[see ([29]) ] and the q/e 1 asymptotic, the latter straightforwardly obtained from 
([281) : Cos « 2(e/g), iVmax ~ 0.11(g/e)^/(l + r/Z) (but this is purely formal because 
the long-time conditions Stg ^ 1 and kyPiStg 3> 1 will be broken in this regime, 
except for extremely long wavelengths). The discrete points show A^max obtained via 
an AstroGK numerical parameter scan: fc||t;thi/*S' = 0.5, varying k±pi and holding all 
other parameters fixed as in figure |4] (note that the asymptotic results do not in fact 
depend on k\\ or S, although the quality of the long-time asymptotic does). 



where we have used (1251) for to and computed the integral numerically. It is clear from 
( 127]) that the integral converges on its lower limit and is not dominated by it, so it does 
not matter that we cannot technically use ( 126|) for short times. 

According to (!28|) . depends on both wavenumbers via us only (a plot of this 
dependence will be given in figure |8]). Assuming o/e ^ 1, we find that the amplification 
exponent is maximised for cos ~ {^/qY^^, givingj 

for Kp,^(-Y''^. (29) 



1 + t/Z \qj S 

These results are illustrated in the left panel of figure [51 The amplification time for the 
maximally amplified modes identified in f[2^ is, from f[2S[) . 

to - ^ ('^Z^)^^^ (30) 

|fc|||^^thi(l + r/Z)y/7c \SkyPi\{l + t/Z)^' 

where we have restored dimensions to make explicit the dependence of to on ^|| • 

Thus, we have learned that an entire family of modes, characterised by a particular 
(linear) relationship between ky and k\\^ given in ( 1291) . will eventually enjoy the same net 
amplification, even though, as follows from the results of § 13. 1[ they were not the fastest 
initially growing modes and some within this equally amplified family started off growing 
more slowly than others or even decaying. The more slowly growing modes are the 

^ Note that this is well outside of the wavenumber domain populated by the damped sound waves, 
(Ds < t/q (see [Appendix B.l[ ). 
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longer- wavelength ones and, according to ([30]), they compensate for their sluggishness 
with longer growth times. 

Note that ( 130|) confirms that the long-time limit is analytically reasonable because 
for large g/e, ( 130|) formally satisfies kyPiSt^ ~ {q/^Y^^ ^ 1 and also StQ ^ 1, provided 
k\\Vthi/S -C g/e. The latter condition is marginally broken by the fastest initially growing 
modes: indeed, in § 13.11 we saw that they had k^\Vthi/S ~ O.lg/e, so for them, Sto ~ 3, 
not really a large number. For all longer-wavelength modes, to is safely within the 
domain of validity of the long-time limit. 

Note also that, as we show at the end of Appendix B.l , the time-dependent 
dispersion relation (|26l) and its consequences derived above can be obtained directly 
from the PVG-instability dispersion relation f lT6|l simply by restoring its dependence, 
setting kx = Skyt, taking the short-wavelength and long-time limit {k^Pi ^ 1, S't ^ 1), 
and assuming e/g ^ <^ 1. This calculation underscores the fundamental simphcity 
of the physics of the transient amplification: perturbations initially destabilised by the 
PVG are eventually swept by the perpendicular velocity shear into a stable region of 
the wavenumber space. 



3.2.3. Limits on short and long wavelengths. We have seen that modes with parallel 
wavenumbers up to A;||ftht/5' ~ g/e can be transiently amplified. From ( l29l) . we conclude 
that of these, the maximally amplified ones will have perpendicular wavenumbers up to 
kyPi < (g/e)^^^, i.e., kypi can be relatively large — unlike in the short-time limit treated 
in § 13. H where the modes with kyPi ~ 1 grew the fastest (although large kypi were also 
unstable). 

It should be understood that, while there is no ultraviolet cutoff in our theory 
that would limit the wavenumbers of the growing modes (in either direction), such a 
cutoff does of course exist in any real system. In the parallel direction, those k\\ that 
were strongly damped in the short-time limit (see flB.l|) ) are unlikely to recover in the 
long-time limit. In the perpendicular direction, the cutoff in kyPi will come from the 
collisional damping, which, in gyrokinetics, contains a spatial diffusion (see, e.g., P), 
and from the electron Landau damping, which we have lost by using the Boltzmann 
electron response (see ([3])) and which should wipe out large ky and k\\. 

On the infrared (long-wavelength) side of the spectrum, we have no cutoffs either. 
In a slab, these would be provided by the dimensions of the periodic box. In a real 
plasma, the cutoffs are set by the scales at which the system can no longer be considered 
homogeneous (in a tokamak, these are the equilibrium-gradient scale lengths and the 
minor radius for the perpendicular scales and the connection length qR for the parallel 
scales; we will need these considerations to fix transport scaling in §[6]). 

3.2.4. Significant amplification threshold. If we maximise fl28l) without assuming 
g/e ^ 1 (with the caveat that the long-time limit asymptotics are at best marginally 
valid then), we obtain a more general curve than (!29|) . plotted in the right panel of 
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figure El We may define a critical threshold for significant amplification: A^max = 1 
when g/e ^ 7. The role of this threshold will be discussed in § 16.11 



3.2.5. Long-time decay. Finally, we obtain the long-time asymptotic decay law. Let 
us seek a solution of ( 126|) such that t 3> to and 7^—1. Then 



-^20F|7|e^' ^ 7^ -Win-. (31) 
to V ^0 

Since 7 is only root-logarithmically large, the quality of this asymptotic is rather poor. 

If we insist on a more precise decay law, we can retain small corrections in fl3T|) and 

get what turns out, upon a numerical test, to be a reasonably good approximation (see 

figure HI) : 



Lt^^t^^ (32) 

For the maximally amplified modes (see fl29|) ). the dimensional damping rate is 7(t) ~ 
I ^11 \vthi'l{t) with to given by ( l30l) . This tells us is that the decay is just slightly faster than 
exponential at the rate of order S. The longest- wavelength modes decay the slowest, 
after having being amplified the longest. 



4. Solution including ITG 

Let us now generalise the results obtained in §|3]to include non-zero (i.e., non- negligible) 
density and temperature gradients. This means that we restore the terms involving 
and Tji in the general integral equation f lTT]) . 



Short-time limit: the ITG-PVG dispersion relation 

The short-time limit, introduced at the beginning of § 13.11 for the case of pure PVG, 
is treated in an analogous fashion for the general ITG-PVG case. An analysis of the 
solutions of the resulting dispersion relation is useful in that its results assist physical 
intuition in ways relevant for some of the forthcoming discussion, but it is not strictly 
necessary for us to have them in order to work out how transient growth happens in the 



presence of the ITG drive. We have therefore relegated this analysis to [Appendix B.2 



^.2. Long-time limit 

We now continue in the same vein as in § 13.21 and consider the long-time limit {St ^ 1, 
kyPiSt ^ 1), in which we can simplify the kernels involving the Bessel functions in (fTTll 
by using (IT^ and also 

A(A.A0.1Wv^-v^)M_gl^. (33) 

^ ' ^ 2 2 2 4 ^ ^ 
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This allows us to rewrite f llip in the form that generalises (I20p : 



1 (cu| + l)A^2^^ 



^{t-M). (34) 



2 4 

As in § 13.21 we seek a solutions to this equation in the form taking At t and 
expanding the delayed potential under the integral according to fl22l) . The result is the 
generalised form of fl23l) : 



2W* 



(35) 



Using again the plasma dispersion function f[T5|) to express the time integrals and 

it) = ^l{t)/ 

- riiU^u(t) 



introducing the complex scaled frequency ui{t) = i'y{t)/ y^a;! + 1, we get 

{q/e)us - 1 



1 + -) V7r((I;| + l)|a;5|t 



[l + u{t)Z{u{t))] 



+ (r/,-l)a;,Z(£l;(t)), (36) 
a time-dependent dispersion relation, which is the generalisation of ( l24l) . 

4-2.1. Transient growth. The general argument that the real part of 7(t) (i.e., the 
effective time-dependent growth rate) must eventually decrease and so fluctuations will, 
in the end, decay, applies to ( 135|) similarly to the way it did to ( l23l) (see § l3.2.ip . although 
this decay need not (and, as we will see, will not) be monotonic. The time to when the 
transient growth ends is now determined as follows. Let Ima)(to) = and a)(to) = ^^o 
(real!). Then from ( 136|) taken at t = to, we find the real frequency ujq by demanding that 
the imaginary part of the right-hand side vanish — this means that the coefficient in 
front of Z{uq) must be zero, because, for real Uo, ImZ(d)o) = y/7ie~'^o. This condition 
gives 

^■wo - (?7i - l)w* = 0, (37) 



- ~2 il/^Ps 



whence 



Uo 



iq/e)us - 1 ± ^\iqle)GJs - \f + 4r/,(r7, - 1) GjI (a;| + 1) 



(38) 



2?7ju;* a/w| + 1 

Substituting this solution into the real part of f l36|) and taking advantage of the already 
enforced vanishing of the coefficient in front of 2(a)o), we get 

, ^ (g/6)a;5 - 1 + ^\{qle)us - 1]^ + r/^l - l/r].) cD^ ((^|TT) 
° 2(l + r/Z)v^(a;| + l)|u;5| ' ^ ^ 

where we have replaced r/j^o)^ = rfgu'^g/A with 775 = Vthi/LxS a new parameter that 
measures the strength of the ITG drive relative to the velocity shear. Note that we 
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picked the "— " mode in ( 155]) because the "+" mode is not amphfied (to < 0, assuming 
rji > 1). Equation f p9|) is the generahsation of f l25|) . to which it manifestly reduces when 
T]s = and with which it shares the property that the transient growth time depends 
on ky and only via us and the time normalisation factor |/i;|||fthj- 

4-2.2. General dispersion relation. We can now recast the general ITG-PVG case in a 
form that shows explicitly how it reduces to the case of pure PVG drive studied in § 13.21 
First we note that the transient growth termination time f p9|) can be rewritten as 



(PVG) _ ^5^5V^5+ 1 

{q/e)ujs - 1 



where 4^^^^ is given by (ES]), rjsujs = 2?7iW*, r]s = Vtu/LrS and cr = sgn[{q/e)us - 
Then the time-dependent dispersion relation (l36l) can be manipulated into the following 
form: 




f = (2 - x^)[l + CuZiu)] + f 1 - - ) X2i^)- (41) 



The analogous equation for the case of pure PVG drive, (1261) . is recovered when x <^ 1, 
which means rjs ^ q/e (cf. f fT3|) ) and r75a;5 ^ 1. In this limit, the behaviour of 
fluctuations in the presence of both PVG and ITG drives is well described by the results 
of § 13. 2[ Before discussing the general case, it is useful to consider the opposite extreme 
of weak velocity shear. 

4.3. Case of weak shear 

Let r]s ^ q/e and r]sOJs ^ 1, so x ^ 1 (note that the same limit is also achieved for 
Us — 00). Then the x dependence falls out of PT]) : 

— = -Co[1+u)Z{Cj)\ + Z{Cj), (42) 
^0 

where we have discarded the l/r/j terms by assuming r/j ^ 1. The transient growth 
termination time in this limit is, from ( HU]) . 

to = - X tr^^ = (43) 

2^ ' 2(l + r/Z)v/vr(a;^ + l) ^ ' 

When t ^ to, we find the solution of ( l42|) by expanding in a; 1. It turns out that it 
consists of a large real frequency and an exponentially small growth rate: (H2|) becomes 

to 2ci) 2t \y1t J 

More generally, for finite values of t/to, the solution of f H2l) is cD = u(t/to), with a 
functional form independent of the parameters of the problem. This solution it plotted 
in figure El It turns out that at t ^ 0.15 to, the growth rate increases sharply, reaches a 

If {q/e)u)s < 1, tg^^'^'' < and x < 0, but the ITG mode can still have transient growth, tg > 0. 




Figure 6. Time evolution of the effective growth rate and frequency: the red 
(upper bold) line is "f{t) = lmuj(t)/\k\\ |uthi\/w^ + 1 and the blue (lower bold) line 
is ReLLi(t)/|fc|| juthi^/wg + 1, both obtained as a numerical solution of (|42t and plotted 
vs. t/to (the time axis is logarithmicio). The black (thin) lines show the growth rate 
and frequency given by the asymptotics (|44p and (jC.9p (the latter taken in the limit 
X — >■ oo and r/i oo). 



finite maximum and then decreases towards zero, which it reaches at t = to, whereupon 
growth turns to decayEI 

Thus, there is a period of strong transient amphfication, which lasts for a finite 
fraction of time to- The amphfication exponent is 

N= dm) = tJul + l [ d^lmuiO - 0.057 (45) 
Jo ^ Jo l + r/Z 

where we have used ( H3l) and calculated the value of the integral under the curve in 
figure [6] numerically (note that since the growth rate is exponentially small at t to, 
the precise lower integration limit is irrelevant). Remarkably, unlike in the case of the 
PVG drive (see fl29|) ). the amplification exponent has no wavenumber dependence at 
all. Also unlike in the PVG case, it does depend on the shear and on the temperature 
gradient: N oc rjs = Vthi/LTS. 

To recapitulate, we have found that, at low velocity shear, all modes are amplified by 
a large (and the same) factor before decaying eventually. Their transient amplification 
time is given by f l43|) . Restoring dimensions, f H3|) and f H5|) are 

W(^r5) AT ^ 0.057 (46) 



2{1 + T/Z)^n{S^k^^p^ + klviy ' 1 + ^/^ 



^ This implies that perturbations first grow due to the ITG-PVG instability (at St -^i 1), then slow 
down to exponentially small growth rates, then (at St ^ 1) grow vigorously again before finally starting 
to decay at t = to- The intermediate period of virtually zero growth, which is a feature both of the 
weak-shear regime and of the general case (see ? l4.4p and may appear strange at first glance, can be 
traced to the kx dependence of ITG-PVG growth rates at long parallel wavelengths — this is explained 
in [Appendix B.2| 



Suhcritical fluctuations in rotating plasmas 



18 




t 
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Figure 7. Effective normalised growth rates lmuj{t/tQ) (red, top) and frequencies 
Keult/to) (blue, bottom) obtained via numerical solution of (j4ip with r]i = 5 and 
X = 0.1, 1, 2, 10 (from top/lighter to bottom/darker curves). See figure [CT] for a more 
detailed depiction of the X = 1 case (in [Appendix C where the functional form of 
these curves is derived analytically). 



The transient growth lasts for a very long time at low 5* and the longest-growing modes 
are the long-wavelength ones. The limit — >^ is singular in the sense that for arbitrarily 
small but non-zero S all modes eventually decay, while for S = 0, the indefinitely growing 
linear ITG instability is recovered (to = oo, N = oo). 

We have already made the point (in § l3.2.3p that while our theory does not limit 
the transiently growing wavenumbers from above, a fuller description of the plasma will. 

4-. 4- Case of finite shear 

In the intermediate regimes between large and small x (i-^-, weak and strong shear), the 
solutions of pT]) transit from the weak-shear form described in § 14.31 to the pure-PVG 
case treated in § 13.21 Figure [7] shows the time-dependent growth rates and frequencies for 
several values of x- As x decreases (i.e., S increases), the peak of the growth rate moves 
further into the past and the growth rate asymptotes to the pure-PVG case (figure H]). 
It is not hard to convince oneself analytically that this is indeed what ought to happen. 
Since intuitively it is rather obvious, further asymptotic considerations on this subject 



are exiled to Appendix C The long-time decay asymptotic is also derived there (it is 
exactly analogous to that found in § l3.2.5p . 

Similarly to our previous calculations, the amplification exponent is 

N{us) = to{cos)^oul+l d^lmu{^,x{cos)), (47) 

where Co is the solution of fj4T]) . The wavenumber dependence enters via the us 
dependence of to and of x (see fj40|) ). The numerically computed amplification exponent 
as a function of us and of rjs is plotted in the left panel of figureOfor r]i = 5 and q/e = 10 
(these are representative of the values encountered in the more realistic numerical studies 
of tokamak transport [3 |25l EE]). The pure-PVG case treated in § 13.21 remains a good 
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Figure 8. The amplification exponent TV, given by (|T7)) with uj the solution of (|^T|) for 
7^,; = 5, g/e = 10 and t/Z = 1. Le/t paneZ; N vs. for 77s = 1, 2, 5, 10, 20, 50, 100, 200 
(from bottom to top curves, in darkening shades of red). The numerical results 
for the same cases are shown in figure [9l The rjs — case (pure PVG drive) 
was virtually indistinguishable from 775 = 1 when plotted (not shown here). The 
lowest-Tys behaviour is well described by ([28|). the highest-775 by (j45|) (in the latter 
case, except for corrections associated with finite 7?^, which are easy to compute if 
they are required). Right panel: The maximal amplification exponent N vs. the 
normalised shear {rjg^ — SLr/vthi)- The maximum A^(ws) is reached at (I15 « 0.54, 
independently of rjs- The dotted lines are the asymptotics ([29|) and (|45|). The finite 
offsets between the asymptotics and the exact curve are due to the fact that the 
asymptotics were calculated in the limits rji ^ 1 and q/e 3> 1, while for the exact 
solution we used relatively moderate values of these parameters. The discrete points 
show iVmax obtained via an AstroGK numerical parameter scan varying kj_pi and S 
(i.e., rjs) while holding k^^Lx = 0.02 and the rest of the parameters fixed at the same 
values as quoted above. 



approximation up to values of rjs of order 10. After that, there is a transition towards 
the weak-shear hmit (§ l4.3p . accompanied by the loss of wavenumber dependence as rjs 
is increased to values of order 100. The constant of proportionality between ky and 
k\\ for the maximally amplified modes (i.e., 0)5 at which is maximised) does not 
appear to depend on rjs, although at large r]s, the maximum is increasingly weak. The 
maximal amplification exponent is plotted in the right panel of figure [HI It is perhaps 
worth pointing out the qualitative similarity between this plot and figure 1 of Ref . [25] , 
obtained from gyrokinetic simulations in full tokamak geometry. 

Finally, the amplification exponent as a function of ky and k\\, obtained in direct 
(linear) numerical simulations, is shown in figure [9] (the parameters are the same as 
in the "theoretical" figure El left panel). The transition from the PVG curve fl29|) 
to the fiat wavenumber dependence ( l45l) is manifest, as are the limits of applicability 
of our approximations in the wavenumber space. Note the different normalisation of 
the parallel wavenumber here (fcyi^T, characteristic of ITG) compared to the left panel 
of figure |5] {k^\Vthi/S, characteristic of PVG). Hence the drift towards higher k\\LT as 
Vs = Vthi/SLx decreases towards the PVG-dominated regime, where the parallel scale 
of maximally amplified modes is set by the shear rather than the temperature gradient. 




Figure 9. The amplification exponent N vs. k\\LT and kypi^ obtained numerically 
using AstroGK for the same parameters as the curves in the left panel of figure [51 
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5. Qualitative summary of the linear results 

In a gyrokinetic plasma with radial gradients of temperature and parallel velocity, both 
gradients are sources of free energy and so will drive the growth of fluctuations (ITG 
and PVG instabilities). The typical growth rate is of order 7 ~ fthi/i^r for ITG and 
7 ~ qS/e for PVG (see ( |T8|) ). or the mean square of the two if they are comparable (see 
(IB.lOp ). Because the mean plasma velocity is toroidal, it always has both a parallel and 
a perpendicular component (the latter a factor of g/e smaller than the former). The 
shear in the perpendicular {E x B) velocity is stabilising and causes all modes to decay 
eventually, so the fluctuation growth is transient — it is always transient in the limit, 
considered here, of zero magnetic shear and it is transient for large enough velocity shear 
S when the magnetic shear is flnite [351 IZ]- If the linear physics provides sufliciently 
vigorous and lasting ampliflcation of finite initial perturbations, it is intuitively clear 
that the system must be able to sustain nonlinearly a saturated (subcritical) turbulent 
state (see §[6]). Therefore, the interesting question is how much transient amplification 
should be expected to occur and on what time scale. 

In the preceding sections, we have addressed this question mathematically, with the 
results summarised by figures O [HI and [9] (see also (|29l) . ( |30l) and (H6i) ). Very roughly, 
these results can be explained as follows. The effect of the perpendicular shear is to 
produce a secular increase with time of the radial wavenumber, kx{t) ~ Skyt. When 
this becomes large enough, the instability is killed by Landau damping (see discussion 



at the end of Appendix B.l ). If we estimate that this happens after to ~ "S*^^ (i.e., for 
kxito)pi ~ 1, assuming kyPi ~ 1), we may conclude that initial perturbations will be 
amplified by a factor of e^, where the amplification exponent is 

iV ~ 7^0 ~ for ITG and iV ~ - for PVG. (48) 

LtS e 

Thus, the shear quenches the ITG amplification — but cannot fall below the shear- 
independent level associated with the PVG (figure [8], right panel). This is indeed the 
case (see (12^ and (H^ ). although, strictly speaking, one has to take into account the 
dependence of the quenching effect on the perpendicular and parallel wavenumbers — 
long- wavelength modes grow more slowly, but for a longer time; in the case of PVG, there 
is also a preferred relationship kyPi ~ {e/qY^^k\\Vthi/ S for the most strongly amplified 
modes (see § 13.21 and figure [9]). While these wavenumber dependences are likely to be 
important in the analysis of the resulting turbulent state and the associated transport, 
they effectively cancel out in the expression for the amplification exponent (because 
7 oc kyPi, to oc l/kyPi at long wavelengths) and the results of the qualitative argument 
that we have given hold true. 

It is instructive to compare these results with the conclusions of a long- wavelength 
fiuid ITG-PVG theory presented in [35] (for the case of finite magnetic shear). In 
that regime, perpendicular shear, by effectively increasing k^it), also caused eventual 
damping of the fiuctuations, but this time via collisional viscosity. Therefore, to estimate 
the transient growth time to, one must set 7 ~ z/jjA;^.(to)p^ ~ UaS'^kyp'^tl, where Ua is 
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the ion collision rate. Then, ignoring wavenumber dependences again, to 7^/^5' ^, so 
the amplification exponent is 

^3/2 I 

N -fto oc -— (X - for ITG and oc for PVG. (49) 

Thus, the E x B velocity shear again quenches the ITG instability, but once S is large 
enough for the PVG drive to take over, the amplification exponent actually grows as 
a/S', the result obtained rigorously by |35] — in contrast with the shear-independent 
N q/t that we have found in the kinetic regime. The practical conclusion from this 
is that it should be easier to obtain states of reduced transport [H [25], [26l [37] in weakly 
collisional, kinetic plasmas. 



6. Subcritical PVG turbulence and transport 

While a detailed nonlinear theory is not the primary objective of this paper, we would 
like to give a preliminary, very crude and very heuristic discussion, inspired by the linear 
results presented above and by previous work on turbulence and transport scalings [8]. 



6.1. Turbulence threshold 

Consider a situation when initial perturbations can grow transiently. The reason they 
decay eventually is that their radial wavenumber kx{t) gets swept by the E x B shear 



from the unstable to damped region (see discussion at the end of Appendix B.l and at 



the end of [Appendix B.2] ). If nonlinear interactions can scatter the energy from these 
modes back into the unstable region before they decay to small amplitudes beyond the 
reach of nonlinearity, then they can be transiently amplified once again and so on. Thus 
a nonlinear saturated state can be sustained — the subcritical PVG turbulence. This 
argument is entirely analogous to the standard paradigm for subcritical turbulence in 
hydrodynamic shear fiows 

The typical time scale for the nonlinear interactions is the nonlinear decorrelation 
time ~ l/k±6vE, where 6ve ~ k±{ap/B) is the fiuctuating E x B velocity. Therefore, 
in order to sustain turbulence the transient growth should last at least as long as one 
nonlinear decorrelation time: 

to > TT^, (50) 
k±dvE 

where to is the amplification time. In the saturated state, if it is sustained, the rate of 
amplification should be comparable to nonlinear decorrelation rate: 

A^ 

7eff ~ T k±6vE, (51) 

^0 

where A^ is the amplification exponent. Combined with (!50|) . this immediately implies, 
unsurprisingly, that the criterion for onset of turbulence is 

N>1. (52) 
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For subcritical turbulence, this criterion replaces the marginal stability condition 
7 = 0, usually employed for cases with well-defined eigenmodes. It is the "significant 
amplification threshold" introduced in § I3.2.4[ As we saw there, for PVG turbulence 
in a slab, it is equivalent to q/e > 7. One might expect that much more sophisticated 
criteria could be derived by refining our arguments and testing these refinements against 
dedicated numerical parameter scans. The key conclusion is that q/e is now the critical 
parameter to be tuned and so magnetic configurations in which it is smaller may hold 
the promise of reduced or even completely suppressed ion turbulent transport (e.g., 
spherical tokamaks [39]). 

6.2. Transport scalings 

The standard mixing-length heuristics (often, somewhat misleadingly, referred to as 
"quasilinear" theory) are based on the argument that if fiuctuations are driven by a linear 
instability with a characteristic growth rate 7, then they will saturate at amplitudes and 
scales such that the nonlinear decorrelation rate is comparable to this growth rate, viz., 
k±SvE ~ 7- Then the turbulent heat diffusivity is xt ~ SvE/k± ~ and so the ion 
heat fiux is Qi ~ xr^iTi/ Lt. 

For subcritical fiuctuations, there is no definite 7, but it is intuitive to argue, as we 
did in §[6], that it should be replaced by 7efr ~ N/to. For the PVG-driven fiuctuations, 
we showed in § 13.2.21 that, at maximal amplification, N q/e and to ~ {q/e)/^Vthi, so 
(I5T1) gives 



This simple estimate is actually consistent with a very general idea that in systems with 
parallel propagation (or particle streaming) and perpendicular nonlinearity, turbulence 
tunes itself into a critically balanced state, viz., the time scales for these two effects are 
always comparable [211 [121 HH HSl [Ml [8] . We may now use the relationship ( l29l) between 
and ky for the maximally amplified modes to estimate k±pi ~ {^/qy^^k^^Vthi/ S and 
conclude, therefore, that 



In the last expression, we have made another important assumption: since it is the 
longest- wavelength fiuctuations that dominate transport, we should use the lowest 
parallel wavenumber possible in a tokamak: k\\ ~ ^/qR, where R is the major radius. 
This prescription was proposed in [8] for ITG turbulence. In a sense, it is even more 
natural here than it was there because, in the theory developed in the preceding 
sections, the wavelengths of maximally amplified modes are not limited from above 
by any microscale physics — if we took the slab model literally, the limit would be the 
periodicity length of the box; in a tokamak, the connection length is a natural choice. 



7eff ~ ^ll^^thi ~ kj_6vE. 



(53) 




(54) 
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Finally, under this scheme, the ion heat flux scales af 



Q^ XT fSVq'/^R 

I I 



(55) 



One fairly obvious feature of the scalings ( 15^ and ( l55l) (independent of most of 
the specific assumptions that we made in deriving them) is that the heat diffusivity 
is independent of the temperature gradient and so heat transport is not very "stiff" 
— the relevant comparison is with the scaling for the ITG regime, Qi oc q{R/LT)^ 
[8] . A softening of transport in the presence of velocity shear has indeed been reported 
both in experimental [321 [33] and in numerical [TJ |25l |26] studies. Physically, it is not 
surprising: as the driver of the turbulence in this regime is the PVG, not the ITG, 
steeper temperature gradients do not produce stronger turbulence and so the positive 
feedback loop between R/Lt and the heat diffusivity is broken. 
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^ There is a number of reasons to take these specific scaling predictions with a grain of salt. Besides 
making the assumptions stated above, we have ignored many effects that may be important and may 
change our estimates of the relevant scales, times and amplitudes: the role of zonal flows in regulating 
and/or sustaining the turbulence, the possibility that the effective radial and poloidal wavenumbers, 
and ky, are not the same in a system with imposed velocity shear, the role of magnetic shear if 
it is present, various geometry (curvature) effects etc. It is also not necessary, although intuitive and 
possibly supported by numerical evidence [26], that the relationship ((29)) between fc|| and ky obtained 
by maximising the linear amplification should persist in the nonlinear regime, especially if the system is 
far above the significant amplification threshold (ii l6.ip . There are several possible alternative theories 
that can be constructed in a similar vein to that presented above, but the current state of numerical 
and experimental evidence does not yet allow us to differentiate between them in a falsifiable fashion. 
Future investigations will focus on this task. Preliminary numerical studies suggest that far from the 
significant amplification threshold f§ l6.1|) , the ion heat flux does not in fact scale as strongly with q as 
suggested by ([55]) (E G Highcock 2011, unpublished). 
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Appendix A. Low-Mach-number local gyrokinetics in a rotating 
axisymmetric plasma 

Here we describe the version of the gyrokinetic system of equations appropriate for a 
rotating axisymmetric plasma, which is the starting point for our calculation. For a 
detailed derivation, we refer the reader to p] (earlier treatments are [H SSj 138]). 

We consider the axisymmetric rotating equilibrium and work in the subsidiary low- 
Mach-number limit, as described at the beginning of §[2l The distribution function of 
particles of species s is written in the following form 
Zseip{r) 



fs{r,v) 



1 - 



Fos{'ip{Rs),£s) + Fu{Rs,es,iJ,s,o-\\) + hs{Rs,es, iJ,s,<y\\), (A.l) 



Let us explain the numerous notation that appears here. The standard 6D kinetic 
position-and-velocity phase space (r, v) is transformed to the 5D gyrokinetic phase 
space, where the dynamics are averaged over the Larmor orbits and so do not depend 
on the gyroangle. If w = v — u is peculiar velocity with respect to the mean flow, B is 
the magnitude of the mean magnetic field, b = B/B its direction, Qg = ZgeB/msC the 
cyclotron frequency, ZgC is particle charge {Z^ = —1), nig particle mass, c the speed of 
light, then the gyrokinetic variables {Rs,es, fis,o'\\) are defined as follows: the guiding 
centre position R^ = r — bxw/Qs, the energy variable Eg = mgw'^ /2 (this is only correct 
in the low-M limit), the magnetic moment = 'msw'j_/2B, and the sign of the parallel 
velocity = w\\/\w\\\ (the subscripts _L and || refer to the mean-field direction b). In 
flA.ip . the particle distribution function is split into the mean Maxwellian, which can 
be shown to depend only on the flux label ip{Rs) via the mean density Ug and mean 
temperature T^, namelyjil 

-1 3/2 



FM{Rs),es) = nsi^iRs)) 



2ttTMRs)) 



exp 



TsiHRs)) 



(A2) 



the mean perturbed distribution function Fig, which contains coUisional (classical and 
neoclassical) effects and will not concern us here, the Boltzmann response containing 
the perturbed scalar potential ip, and the guiding-centre distribution function hg. 

If we take the mean fields B (i.e., ip and /(V^)), ns{ip), Ts{ip) and (^{ip) to be known 
and the fluctuations about them to be purely electrostatic, then the latter are fully 
described by a closed system containing the gyrokinetic evolution equation for hg and 
the quasineutrality condition determining 



dhs 
dt 



+ [u{R,) +wiib + Vd+ {Ve)r, 



dR. 



ZseF 



d 



-((Ve)k. ■ VV^) 



dlnng 

dip 



Es 3 

+ ' t; ~2 



dip 



+ 



dR, 

B^ nisWwRduj 
'b'^Ts # 



Fos, (A.3) 



^ Establishing the Maxwelhan equilibrium depends on the plasma being sufficiently coUisional, namely, 
that the collision frequency is not smaller than the fluctuation frequency by more than one order in 
the gyrokinetic expansion parameter. The density being a function solely of ^ only holds in the low-M 
limit. 
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E ^1?^) V = J2^^'I d'^(h,)r. (A.4) 

where w\\ = (T|| A/2(es — fJ-sB), = I{ip)/R is the azimuthal magnetic field, Vd = 
{c/ZseB)b X (msWj|6 ■ V6 + ^is^B) is the magnetic drift velocity, Ve = {c/B)h x Vy? 
is the E X B velocity, and the gyroaverages are defined 

1 Z'^'' 

{v)rs = — / ddif{Rs + 6 X w/Vts) (function of Rs), (A. 5) 

2vr Jo 

1 f'^'" 

{hs)r = — / di!} hs{r — b X w/Qs,£s, fJ'syC'w) (function of r and lu). (A. 6) 

27r Jq 

Note that we suppressed the collision term in (1A.3I) (formally, we are ordering the 
collision frequency small via a subsidiary expansion). 

The approximation of Boltzmann electrons amounts to setting /ig = in f lA.4p . 



Appendix A.l. Local Cartesian frame in the case of zero magnetic shear 

Let us take the mean magnetic field to be locally straight and uniform, so it has constant 
magnitude, no curvature and no shear. This means that the magnetic drifts vanish 
{Vd = 0) and we can introduce a local orthogonal Cartesian coordinate system (see 
figure EI) : 

JDprL JDpIt 

where Bp = I'Vipl/R is the poloidal component of B. We will view this coordinate 
system as having its origin {x = 0) at some reference fiux surface ipo, so in the vicinity 
of this fiux surface, ip ^ ipo + xBpR. Then 



ujR + xBpR^^ 

dip 



§i + §^),(A.8) 



B B 

where all quantities in the last expression are taken aX ip = ipQ and B^ = I{ip)/R is 
the azimuthal component of B. We have only retained the spatial dependence of u (in 
the form of a Taylor expansion) because we are formally ordering velocity gradients as 
0(1/M), while all other mean fields are assumed to have 0(1) variation on the system 
scale. The constant part of the velocity can be removed by going to a frame moving 
(azimuthally) at this velocity. The z component of the velocity shear can be neglected 
because it multiplies d/dz in flA.SP and we order x as small (in the gyrokinetic expansion 
parameter) compared to the parallel scale of the fiuctuations, whereas the y component 
of the shear remains important because it multiplies d/dy ^ d/dz. It follows that the 
velocity field can be replaced by the linear shear fiow as stated in ([1]). 

Given the shearing box approximation described above and the assumption of 
Boltzmann electrons {h^ = 0), the conversion of ( ]A.3I) and ( ]A.4I) into ([2]) and ([3]) is 
straightforward. Note that {c/B)(Ti/Ze) = v^Y^j2Qi = piVtu/l, where pi = Vtu/^i is 
the ion Larmor radius. 
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Appendix B. PVG and ITG-PVG dispersion relations 

Here we provide some elementary analytical considerations of the linear dispersion 
relations that govern the early evolution of the fluctuations. 



Appendix B.l. PVG dispersion relation 

In the pure PVG case, the dispersion relation is ( |T6|) . containing the well-known PVG 
instabihty [11] and the sound wave. The following considerations provide some analytical 
support for figure [3l 



Marginal stability threshold. This is easy to obtain analytically without any 
approximation: at the threshold, Imo) = and then, from f|T6|) . the real frequency 
is readily shown to vanish as well, so we get {q/e)ijJs = (1 + 't/Z)/To{\). Thus, we have 
an instability if|§ 

fcpth» kyPiVoiX) ^ 0.66 

qS/e 1 + t/Z - 1 + r/Z' ^ ' ' 

where A = kypf/2 (see figure |3|, right panel). The second inequality in (IB.1|) implies 
that there is an absolute finite limit on the parallel wavenumbers at which the PVG 
instability can survive (reached for kyPi ^ 1.26). 

PVG modes and sound waves. The growing modes have no real frequency. In fact, there 
is no real frequency as long as {q/e)u}s > I5 which means that the mode is purely decaying 
between {q/e)u!s = 1 and the marginal stability boundary {q/e)u!s = (l+''"/^)/ro(A) > 1 
(see flB.ip ). For {q/e)u}s < 1? the mode turns into a damped sound wave (this is shown 
in the left panel of figure E]). Thus, the PVG drive has pushed the sound waves into a 
wedge in wavenumber space, 

kyP. < ^ (B.2) 

qS/e 

(see figure El right panel) and populated the rest with nonpropagating growing or 
decaying modes. 



Growth at long parallel wavelengths. Looking for unstable solutions in the long- 
wavelength limit {k\\Vthi/S ^ 1), we consider the asymptotic form of (fT6i) with u = i'j 
(pure growth) and 7^1, when 1 + i'jZ{i'y) ^ 1/27^. Then the growing solution is 



l^^gj)^s^]__ {qS/t)kyp,k\\v,u 
^ ^/2[(l + r/Z)/ro(A)-l] ^ V2[(l + ^/^)/ro(A)-l]- ^'^ 

^ If qS/e > 0, the unstable mode has either ky > 0, fc|| > or ky < 0, fc|| < (see figure IBll and 
discussion in [Appendix B.2| ; if qS/e < 0, then ky and fc|| have to have opposite signs. We can assume 
without loss of generality that all these quantities are positive. 
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In the second expression, we have restored dimensions and assumed {q/e)uJs ^ 1 (which 
ensures 7^1 and is consistent with the short-time ordering adopted at the beginning 
of § 13. ip . Note that, while this asymptotic has a peak at kyPi ~ 1, it does not capture 
the maximum growth rate because of its monotonic increase with k\\. The maximum 
growth rate is, in fact, reached for 7 ~ 1, where the plasma dispersion function does 
not yield itself to a simple asymptotic expansion. The numerical solution is shown in 
figure [31 



Growth at short perpendicular wavelengths. Finally, we note that one can obtain a 
good approximate preview of the results of § 13.21 if one restores the dependence on by 
setting A = {k^ + ky)/2 in (fT6|) . The effect of the perpendicular shear is to increase 
the instantaneous kx of the mode, so the limit of short perpendicular wavelengths 
(or, more precisely, large k^ but finite ky) is similar to the limit of long times. Since 
ro(A) ~ I/VQjtX for A 1, (lT6l) becomes (assuming also kyPi k\\Vtu/ {qS/e), i.e., we 
are far outside the domain of the damped sound waves discussed above) : 

where to is given by fl25l) . but with dimensions restored (i.e., not normalised by A;||t'thi) 
and e/g <^ 6)5 ^ 1 (as will be the case for the most strongly amplified modes in § l3.2.2p . 
This formula has two interesting consequences. 

Firstly, setting k^ = Skyt and S't <^ 1, we recover the time- dependent dispersion 
relation (|26|) and so the transient growth and eventual decay of the fluctuations derived 
more formally in § 13.21 (without the assumption e/q ^ 0)5 ^ 1). Thus, there is a 
smooth connection between the short-time and long-time behaviour of the PVG-driven 
fluctuations. Note that applying the same asymptotics to (lB.3p . we recover (127]) . That 
the PVG growth rate must be extinguished at large enough kx is also obvious from the 
instability criterion fIB.ip . where increasing A eventually breaks the first inequality. 

Secondly, in the short-time but also short-perpendicular- wavelength limit {kx <^ky 
but kyPi 3> 1), the growth rate of the PVG instability is independent of ky, as is indeed 
manifest in figure |3l 



Appendix B.2. ITG-PVG dispersion relation 



In the short-time limit, we can, analogously to the derivation in § 13. ![ reduce the 
integral equation ( fTTj) to a dispersion relation for the normalised complex frequency 

u = uj/\kii\vthi- 



r/Z 



- 1 



ro(A) 

where ro(A) 



■Us 



1 — r]iOJ^:iSj [1 + ujZ{u)] + r/jW* 



(B.5) 



^/o(A) and A(A) 



.2 r]i 
1 — A + A/i(A)//o(A); we recapitulate the definitions 



of A = kyPi/2, LUs = Skypi/k\\Vtu, ^* = kyPi/2\k\\\Ln and r/j = Ln/Lr. This is the 
standard slab ITG-PVG dispersion relation for an unsheared gyrokinetic plasma. 
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Marginal stability thresholds. The way in which ITG and PVG coexist is easiest 
to understand by examining the marginal stabihty thresholds. Setting Imw = 0, 
w = Rew = Wo and demanding that the imaginary part of ( IB. 51) vanish, we get the 
equation for the real frequency of the mode at marginal stability: 



2 Vi. 



0. 



(B.6) 



Substituting the solution of this equation back into fIB.Sp . we arrive at the marginal 
stability condition: 



2-2 



+ 



-iz 



q 1 + r/Z' 



0. 



(B.7) 



L ro(A) 

Since riiU^ = kyPi/2\k\\\LT and ujs = SkyPi/k\\Vthi, the above equation can be solved for 
the two curves in the (/cy, ky) plane that enclose the unstable region: 

kyp.r^l [3/2 - A(A) - / [(1 + t/Z)/T,{X) - 1] 



k\]L. 



Sv^TT^ \-i±V^ + vl [3/2 - A(A) - I [1 - ro(A)/(l + tIZ)\ 



(B.8) 



where L^t = Vfti/ ^J^qS/ef + {vtu/LTY = {qS / e) y^YTr^l is a convenient 
normalisation of the parallel wavenumber for the mixed ITG-PVG regime and rj^ = 
{vthi/LT)/{qS/e) = ris/{q/e) is a parameter that measures the relative strength of the 
ITG and PVG drives. 

The pure PVG regime is r/^ <^ 1, in which case the "+" curve above turns into 
fIB.ip (which is perhaps easier to infer directly from flB.7l) ). while the "— " curve is simply 
k\\ = — these demarcate two symmetric PVG-unstable regions aX k\i > 0, ky > and 
A;|l < 0, ky < (assuming S > 0). The pure ITG regime is rju ^ 1, in which case 
Lt and 



iky Pi 



3/2 - A(A) - 1/r], 



4[(l + r/Z)/ro(A)-l/2r 



(B.9) 



These curves enclose four symmetric ITG-unstable regions in the four quadrants of the 
(/cy, ky) plane. In the general case of finite rju, when both ITG and PVG drives play a 
role, these drives combine constructively in the PVG-unstable quadrants /cy > 0, ky > 
and < 0, ky < 0, whereas in the remaining two quadrants, the stable PVG mode, 
which is just a PVG-modified sound wave (cf. (fT7|l ). has a stabilising influence on the 
ITG drive. 



Solutions. It is not hard to show that the solutions of the ITG-PVG dispersion relation 
(IB. 5 1) can be expressed in the general form, which is a direct generalisation of ( ITSl) : 

w = ^/{qS/ey + {vtu/LrY f [kyPi,k\\LuT\riu,rii) ■ (B.IO) 

The growth rate 7 = Imcu, obtained numerically for a typical situation with finite ?7„, 
is shown in figure IBll as are the marg inal stability thresholds ( iRSll . (iRll and (iRTll . 
Note that in general the ITG-PVG mode also has a real frequency. Note also that the 
solutions of the general ITG-PVG dispersion relation retain the PVG-mode property 
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Figure Bl. The ITG-PVG growth rate -f/^{qS/e)^ + {vtht/Lr)^ vs. fc||L„T and kyp^, 
where L„t = Vthi / y/ {qS / e)^ + (vthi/LrY- In this plot, = {vtY,J Lt) / {qS / e) = 5, 
rji — b and t/Z = 1. Only positive values are plotted, black means 7 < 0. The red 
curves show the stability boundary (jB.8|) . the red dotted curves the pure-PVG (ry^ — > 0) 
stability boundary (jB.lj) and the white dotted curve the pure-ITG (r/„ — )■ 00) stability 
boundary (|B.9I) . See figure [B2] for cuts along the kyPi axis for some representative 
values of k\\LuT- 



of independence of ky at short perpendicular wavelengths. This is easily demonstrated 
analytically by taking the limit kyPi ^ 1, k\\LuT in (IB.SP — all factors of kypi then 
cancel on both sides of the equation. 



Growth at long parallel wavelengths. To assist physical intuition and some of the 
forthcoming discussion, it is useful to consider the ITG-PVG dispersion relation in 
the "fluid" limit, uj ^ 1 (case of long parallel wavelengths). Expanding Zip) = 
— l/oj — l/2(jj'^ — 3/46}^ + . . ., we recast (IB. 51) as a cubic equation 



1 + t/Z 

ro(A) 



u + 2riiU^ 



l-A(A)-- 



A(A) + - 
Vi] 



0. 



(B.ll) 



If we also consider long perpendicular wavelengths, A — )■ 0, and assume r^j ^ 1, we get 



2r 



,-,3 



q - 



1) U — T]iUJ^ = 0. 



(B.12) 



The long-wavelength PVG instability f ll7p is recovered as a balance of the first two terms 
(in the limit <^ 1); the fluid version of the ITG instability obtains from the balance 
of the first and third terms (when rju ^ 1): 



(B.13) 
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Figure B2. Left panel: The ITG-PVG growth rate, calculated numerically from ()B.5[) 
for the same parameters and normalised in the same way as in figure IbTI vs. kyPi. The 
three curves are for k\\LuT = 0.01 (red), 0.02 (brown) and 0.05 (black). Right panel- 
Same, but vs. k^pi. 



(three roots, one real, two complex, of which one unstable). 

Analysing the long-wavelength dispersion relation in great detail is not a useful 
exercise because it turns out to be a very poor quantitative approximation to (IB.Sp in 
most parameter regimes. It does, however, make transparent how the ITG and PVG 
modes coexist. It also helps understand qualitatively what is perhaps a somewhat 
obscure property of the ITG (and ITG-PVG) mode: if the parallel wavenumber is held 
fixed and low {ki\LuT ^ 1), the mode is unstable at long perpendicular wavelengths 
(A ^ 1), becomes stabihsed at finite A, and then reignites at larger A. This is because 
the "fluid" ITG solution f lB.13|) depends on the cancellation A(A — ?■ 0) ~ 1 (and the 
limit rji ^ 1) in the second term of (IB.llI) . As A increases, this second term becomes 
larger than the first and the mode is stabilised (or, more precisely, it still has an 
exponentially small growth rate originating from the iy/n e~'^^ term in Z{cd), neglected 
in the derivation of flB.ll|) — this comes from the Landau pole and is not recoverable in 
fluid theories). The instability is rekindled at larger A because the coefficient in the first 
(cubic) term in ( ]B.11|) grows (~ \/X as A — > oo) and comes back into play (although 
the approximation u ^ 1 breaks down simultaneously). 

This behaviour, calculated from the full dispersion relation ( IB.5I1 . is illustrated in 
figure IbTI In the right panel of this figure, we show the growth rate dependence on k^pi, 
the dependence on which we have restored by setting A = {k1+ky)p1 /2 (similarly to what 
was done at the end of Appendix B.l ). The attenuation and subsequent resumption of 
growth as the perpendicular wavenumber grows is the origin of similar behaviour, but 
as a function of time, seen in the long-time limit of the ITG-PVG transient growth in §|1] 



and [Appendix C| — because the effect of the perpendicular shear is to increase gradually 
the instantaneous kx of the mode. At long parallel wavelengths, as time increases, the 
mode first grows (when St 1), then slows down for a while, then (at St ^ 1) the 
growth is resumed for another period and finally extinguished at t = to- 
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Appendix C. Transient growth with finite shear 

Here we provide some asymptotic considerations that back up the behaviour of the 
time-dependent growth rates and frequencies described in § 14.41 

As in § 14.31 when t ^ to, we may expand the time-dependent dispersion relation 
dH]) assuming large frequency a) ^ 1, which is mostly real (7 = Ima) ^ Rea)): 




1 ^ _ A _ 1") A _ + (2 - ;^cD) (I; e-^^ (C.l) 



The two terms on the right-hand side that involve x ^^e the ITG terms and the other 
two are the PVG terms. When the latter or the former predominate, we recover the 
limiting cases treated in § 14.31 and § 13.21 respectively (in the case of strong PVG-driven 
growth the resonant term should be dropped as the growth rate is not small). The 
transition between the two regimes can be understood essentially by calculating finite-x 
corrections. 

Strong ITG, weak PVG. Let us first consider the case where ^ 1 but x ~ 1- Then 
the two PVG terms on the right-hand side of (IC.ll) can be dropped and we obtain the 
following solution 

^] - '/^-^^ \ 7-^^^^e-^ (C.2) 

l + a^l + (l-l/r^,)x^ 2t' 1-2/r/, 

a straightforward generalisation of fj44|) . to which it reduces when x — )■ 00 and 774 — )■ 00 
(there is, of course, no problem retaining finite rji, it was assumed large in § 14.31 merely 
to make the algebra more compact). The exponentially small growth rate 7 stops being 
exponentially small and takes off to order-unity values when the frequency u stops being 
large and becomes a) ~ 1. It is evident from (lG.2p that this happens at earlier times for 
smaller x- Hence the behaviour of the growth rates and frequencies shown in figure [T) 
the peak growth moves into the past (smaller t/to) as x decreases. The effect of finite x 
is only noticeable when x^ ~ 1; at larger values, the frequency is basically independent 
of X and so the weak-shear limit (x — ?■ C)o; § 14. 3p applies. 

Strong PVG, weak ITG. Let us now restore the PVG terms in (IC.ll) . This becomes 
necessary when x^ ~ 1 or smaller. In this limit and in neglect of the exponentially 
small resonant term, the solution of ( IC.ll) is 



w ~ -4 1 ±i, 



(PVG) 


2t 



(PVG) 



ViJ t 



2 



(C.3) 



where, to make the formula more compact, we have used the first expression in (140 p to 
express the left-hand side of (IC.ip in terms oit'^^'^\ Here we are assuming {q/e)ujs > 1, 
so a = 1, > and X > 0. 
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Figure CI. Effective normalised growth rate Imw(t/to) (red, top) and frequency 
Rea)(t/to) (blue, bottom) obtained via numerical solution of (j4ip with rji = 5 and 
X = 1. This is one of the cases already shown in figure [3 replotted here to compare 
with the analytic solutions derived in [Appendix C[ the black (thin) lines show the 
(combined) growth asymptotics (|C.3|) and (|C.6|) for t/ig ^ (the growth rate (jC.7|) is 
so small that it is zero for all practical purposes) and the decay asymptotic (|C.9p for 

t/lQ > 1. 



The "+" branch in (IC.3I) is a growing mode and is a direct generahsation of ( 1271) . 
The difference with the pure-PVG case is that growth does not extend indefinitely into 
the past, but requires 

2 / 2 



t 



(PVG) 32 

''0 



> ^ I 1 - - 



At longer times, the growth rates increases, reaches its maximum value 

2 / o \ 2 



7n 



at 
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(PVG) 




IQ 



^ 11 - - 



(C.4) 



(C.5) 



and then asymptotes to the pure-PVG case (!27|) . At shorter times, the frequency 
(IC.Sp becomes purely real (asymptoting eventually to u given by (lC.2p ). To get the 
growth rate in this regime, we must restore the resonant term in fIC.ip . which gives an 
exponentially small growth rate. The (growing) solution now is 



7 



^ I 1 - - 



(PVG) 




\ 



^ I 1 - - 



(PVG) 




(PVG) 




2t 



X(l - 2/r7,)4P^^V8tl -to/2t 



(C.6) 
(C.7) 



Figure [CT] shows the numerical solution of f l4T|) together with the asymptotics flC.3|) . 
flC.6|) and flC.7|) . all calculated for x = 1 and r]i = 5. Even though the asymptotics are 
only technically valid for x<l/(I;^l, we see that they in fact work rather well even 
for moderate finite % (this is due to the accumulation of numerically small prefactors 
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that multiply x the above expressions). For values of x significantly larger than unity, 
this agreement breaks down and one has to use flC.2p . 



Thus, we have learned that both PVG- and ITG-dominated dominated transient 
growth is exponentially weak at very short times, then increases to large or finite values 
before slackening again and turning to decay at t = to- At what time the transition 
happens and how large 7 can get depends on x- Recalling the definition of x (see f HOj) ). 
we note that larger values of x ^i^re achieved for stronger shear (i.e., rjs larger compared 
to g/e), larger kyPi or smaller k\\v^^Yii/ S (i.e., larger us). Thus, the transition between the 
ITG and PVG regimes happens non-uniformly in the wavenumber space (as, indeed, is 
evident in the middle panel of figure [8]) . 

We remind the reader that all of the above is only valid in the long-time limit, 
ugt ^ l^kyPi (see the beginning of § 13. 2p . which means that for some wavenumbers 
and/or values of rjs and g/e, the transition from exponentially small to finite growth 
rate may be superseded by the transition from the short- to long-time limit, i.e., the 
initial ITG-PVG instability may not have time to peter out (due to increase in k^pi 
caused by shearing; see Appendix B.2 ) before being rekindled again by the transient 
growth. 



Long-time decay. Finally, to complete our treatment of the general ITG-PVG case, let 
us consider the eventual decay of the fluctuations. This is done in exactly the same 
way as in § 13.2.51 In the limit t 3> to? we assume again that u = ij, where 7 is large, 
negative and mostly real. Then fj4Tl) is approximated by 




t 

To 



2v^(2 + ^x|7l 
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In 



to' 



so, like in f l3T]) . we have a root-log law, i.e., the decay of the modes with time is just 



barely super-exponential. As before, ( 1C.8I) is a quantitatively poor approximation except 
at ridiculously long times and a better one can be obtained by retaining corrections. 
This way we also obtain the (decaying with time) real frequency. The result is 
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In 



l + Vl + (l-l7^ t/(2v/?to) 
v/4 + x2 1n(t/(2v/¥to) ^/Ht/{2^to) 



arctan(x|7l/2) 
2|7| 



(C.9) 



For X ^ 1, we recover the pure-PVG result fl5^ . For x ^ 1 and rji ^ 1 (the weak-shear 
limit), these asymptotics are shown in figure [6] and for x = 1 and r/j = 5 in figure [CTl 
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